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Abstract 

The renormalon calculus is used to calculate the terms of order /Jq -1 ^ in the 
perturbative expansions of the Wilson coefficients and hard-scattering kernels en- 
tering the QCD factorization formula for hadronic i?-meson decays into two light 
pseudoscalar mesons. The asymptotic behavior of the expansions is analyzed, and 
a minimal model of power corrections arising from soft "non-factorizable" gluon 
exchange to the B — * irK, irir decay amplitudes is obtained, which takes into 
account the structure of the leading and subleading infrared renormalon singulari- 
ties. Whereas the resulting power corrections are generally very small, some of the 
strong-interaction phases of the hard-scattering kernels receive sizeable two-loop 
corrections. The implications of these findings on CP asymmetries and branching 
ratios are investigated. 
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1 Introduction 



A theoretical understanding of hadronic B decays is necessary to extract CKM matrix 
elements and CP-violating weak phases from data taken at the B factories. The formal- 
ism of QCD factorization developed in |], |2|, |3| has shown that two-body hadronic B 
decays simplify considerably in the heavy-quark limit. In B — > nK decays, for example, 
QCD factorization states that the matrix elements of the local operators Qi appearing 
in the effective weak Hamiltonian obey the factorization formula 

{irK\Qi\B) = Ff^T^ ® f K $ K + F^ K ® /A 

+ T! 1 ® f B $ B ® f K $ K ® f n $ n , (1) 

where <3?m are leading-twist light-cone distribution amplitudes, and the ^-products imply 
integration over the light-cone momentum fractions of the constituent quarks inside the 
mesons. The hard-scattering kernels Tj are perturbatively calculable quantities. The 
factorization formula is true to all orders in perturbation theory, but it receives power 
corrections in Kqqd/ijii,. If QCD factorization is to be used as a tool for extracting 
the flavor parameters of the Standard Model, it is important to show that these power 
corrections are not anomalously large. Some sources of power corrections were identified 
and estimated in []2], [| |J]. In |5|, ||, the renormalon calculus was used to study power 
corrections due to soft gluon exchange in the decays B — > D^L, where L is a light 
meson. Here we generalize and extend these analyses to focus on soft gluon corrections 
to the amplitudes for B decays into two light pseudoscalar mesons. 

The model of power corrections that we develop relies on the analysis of infrared 
(IR) renormalon singularities, which arise because the low-momentum regions in loop 
diagrams give increasingly important contributions at higher orders in perturbation the- 
ory H || , rendering perturbation theory a divergent asymptotic expansion. IR renor- 
malons are most transparently studied through singularities in the Borel transforms of 
perturbation series with respect to the coupling constant |nj. The Borel transform S(u) 
of the perturbation series for a physical quantity has poles (or branch cuts) at half- 
integer values of the Borel variable u, which lead to ambiguities of order (Aqcb/Q) 2u 
when performing the inverse transform, a process to which we will refer as Borel resum- 
mation. (Here Q is the large momentum scale in the process.) Since the final answer 
for a physical quantity must be unambiguous, these renormalon ambiguities should be 
combined with power corrections of the same order in 1/Q. Therefore, an analysis of 
IR renormalons may reveal non-trivial information about power corrections to various 
processes, which is particularly useful in cases where such corrections cannot be analyzed 



using an operator product expansion (see flfl" , l2| , |13| , [14| , |15| , |16|| for early applications of 
this idea). 

We calculate the Borel transforms of the hard-scattering kernels appearing in the 
QCD factorization formula using the so-called "large-/?o limit" , in which the momentum 
flow through gluon lines is traced by inserting fermion bubbles on the propagators, 
replacing n/ — > — 3/3q/2 at the end of the calculation ||. In this way, the diagrams 
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are effectively computed using a momentum-dependent running coupling at the vertices. 
This procedure gives the dominant higher-order contributions to a perturbation series 
in the very formal limit where n/ — > — oo, in which case /3q — > oo with fixed number 
of colors. Although for real QCD /3 = 11 — |n/ is not a large parameter, the use of 
the large-/9o limit is motivated by the empirical observation that perturbation series are 
often dominated by terms of the form (3q . 

In this paper, we calculate both the Wilson coefficients and the hard-scattering kernels 
entering the QCD factorization formula at next-to-leading order (NLO) in the large-/?o 
limit. This is necessary for obtaining renormalization-group (RG) invariant expressions 
for the physical decay amplitudes. The expansion of the Wilson coefficients in the large- 
(3q limit is a novel feature of our work, which goes beyond the approach taken in |5|, ||. 
The numerical study in Section ^ helps to elucidate the main results of our work. We find 
that the power corrections to QCD factorization due to soft-gluon exchange are small, 
typically of order (Aqcd/^()) 2 - The only exception is a first-order power correction to the 
hard-spectator term, which is to be expected considering that the effective scale for gluon 
exchange with the spectator quark is of order yA^co?^ . Perhaps our most important 
result is that the strong-interaction phases of the decay amplitudes increase significantly 
when going beyond one-loop order in perturbation theory. This can potentially enhance 
the direct CP asymmetries in decays such as B — > tttt and B — > ttK with respect to 
the predictions found at NLO in RG- improved perturbation theory in ||. Whether or 
not such increased asymmetries are realistic depends on whether there are significant 
multi-loop contributions not captured in the large-/?o limit. We leave a more careful 
analysis of such contributions for future work. 

2 Wilson coefficients in the large-/3o limit 

The QCD factorization formula provides model-independent expressions for hadronic 
weak decay amplitudes at leading power in the heavy-quark expansion. All non-trivial 
strong-interaction effects, which give corrections to the so-called "naive" factorization 
(vacuum insertion) approximation, are contained in RG- invariant parameters tij, which 
are products of Wilson coefficients with convolutions of hard-scattering kernels and me- 
son light-cone distribution amplitudes. In order to obtain RG-invariant results, it is 
crucial that the Wilson coefficients and hard-scattering kernels be treated in the same 
perturbative approximation scheme. Since our goal is to resum the perturbation series 
for the hard-scattering kernels in the large-/?o limit, we must first compute the Wil- 
son coefficients appearing in the effective weak Hamiltonian in the same limit. This is 
non-trivial, since the counting scheme underlying the large-/?o limit is different from the 
usual counting of logarithms used in the calculation of the Wilson coefficients. Specifi- 
cally, whereas a s ln(M/ fi) is considered 0(1) in the usual RG counting, this quantity is 
0(1/ Po) in the large-/3 limit. On the contrary, whereas 1/ ln(M//x) is usually treated as 
a small parameter, it is considered 0(1) in the large-/?o limit. 

Consider, as an example, the rare hadronic decays B — > tttt, for which the effective 
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weak Hamiltonian (neglecting small contributions from electroweak penguins) is [HJ 

V^p=«,c ^i=l,2 i=3,...,6 ' 



(2) 



where X p = V p bV* d are products of CKM matrix elements. The scale dependence of the 
Wilson coefficients is governed by the evolution equation (in compact matrix notation) 



C(/i) = T a exp 



da 



a a (M) 



2a /3(a) 



C(M) 



(3) 



where the evolution matrix is given in terms of an ordered matrix exponential involving 
the anomalous dimension matrix 7 and the /3-function. The initial values C(M) at a 
scale M ~ mw are obtained by matching the Standard Model onto a low-energy effective 
theory, in which the heavy fields of the electroweak gauge bosons and the top quark are 
integrated out. 

To evaluate (§) at sub leading order in the large-/3 limit, we expand the QCD j3- 
function as /3(a) = b + 0(1/ Pq), where the rescaled coupling 



b(/i) = fa 



4vr "ln(/iVAy 



(4) 



is 0(1) in the large- Pq limit. Next, we change variables from a to b in the integral in (|3|), 
write the anomalous dimension matrix as a function of b, and expand in powers ofl/Po. 
Taking into account that factors of n/ must be replaced with —3Po/2, we find that in 
2x5 block form the resulting matrix is given by 



l T {b) 












bC 



1 



Mb) 





bB ± 


Ci(b) 



o(i/Pl) 



(5) 



The justification for writing the leading term in this expansion as a constant 5x5 matrix 
Co times the coupling constant b is that this corner of the anomalous dimension matrix 
does not pick up an additional factor of rif when going from leading order (LO) to next- 
to-leading order in perturbation theory, so that all further perturbative contributions 
are absorbed in the matrix Ci(b), the form of which is irrelevant to our discussion. A 
similar argument applies to the submatrix in the lower left corner, which at order 1/Po 
is given by a constant 5x2 matrix B± times b. 

The matrix A±(b) governing the evolution of the current-current operators Q{ 2 in 
the large-/?o limit is a non-trivial function of the coupling, which has been calculated in 
The result is 

A 1 (b)= 1 (b)\ J _M, (6) 
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where N = 3 is the number of colors, and 

26(3 + 26) (1 -6)T(4 + 26) 



7(6) 



(7) 



3r(i-6) r 2 (2 + 6)r(3 + 6) ' 

The matrices Co and B± can be derived from the expressions for the LO anomalous 
dimension matrices given in [H]: 



Co 








1 

N 







1 

N 










-1 







— 1 










1 

N 







1 

A' 










-1 







— 1 





V 


9 


UN , 
6 " t " 


29 


9 


8iV 25 





2 


6iV 


2 


3 6N 



B 1 



2 

~3iV 
2 
3 

2 

~3JV 
2 
3 



o\ 







UN 29 o 
\ 9 9N ° 



Taking into account the textures of these matrices as well as corresponding relations 
between the matching conditions for the various coefficients at the scale \x = mw, we 
find that the solutions for the Wilson coefficients can be expanded as 



c 1 (ri = i + -cl 1 \ f ,) + o(i/ti) J 

Po 



a(/i) = -^cfV) + o(i/A 

Po 



OJ i 



6, 



1 



Po 



where 



-iVCfV) 



cf (m) = C7^(/i) = -^V) = -jvc^o*) 



,(1) 



,(1), 



(9) 



(10) 



Up to 0(1/ Po), the Wilson coefficients C\, . . . , C 6 are thus determined in terms of two 
functions, given by 



6(m) 



( /6 



7W 
26 2 



6(Af) 



f 6M 
U(M) 



Ci 1} (M) 



(11) 



The coefficient of the chromo- magnetic operator Q 8g will be needed only at LO, where 
it is scale independent, i.e., 

Cg ) ( A i) = Gg ) W- (12) 

Note that, as a consequence of the fact that the matrix Co has eigenvalues and 
—2, the scale dependence of the penguin coefficients C3, . . . ,Cq governed by the function 
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Table 1: Comparison of the Wilson coefficients in the large-/3o limit with those ob- 
tained at NLO in conventional, RG-improved perturbation theory. Input parameters 
are m t (m t ) = i67GeV, m b (m b ) = 4.2 GeV, a s (m b ) = 0.2244, and a s (m w ) = 0.f202. 



NLO 


Ct 


c 2 


c 3 


c 4 


c 5 


c e 


fi = m b /2 


1.137 


-0.295 


0.021 


-0.051 


0.010 


-0.065 


/i = m b 


1.081 


-0.190 


0.014 


-0.036 


0.009 


-0.042 


H = 2m b 


1.045 


-0.113 


0.009 


-0.025 


0.007 


-0.027 


Large- /?o limit 


Ct 


c 2 


c 3 




c 5 


c 6 


fi = m b /2 


1.087 


-0.260 


0.020 


-0.061 


0.020 


-0.061 


fi = m b 


1.061 


-0.183 


0.014 


-0.041 


0.014 


-0.041 


H = 2m b 


1.039 


-0.118 


0.009 


-0.028 


0.009 


-0.028 



C^\n) is extremely simple, and given by b(fi)/b(M) — 1 = 26(/x) \n(M/fi). We will 
see below that the appearance of a logarithm is indeed required to cancel the scale 
dependence of the hard-scattering kernels associated with penguin contractions of the 
operators Qi. 

The remaining step in the calculation of the Wilson coefficients is to specify the 
matching conditions at a scale M ~ my/- In general, these will be complicated functions 
of the coupling b(m\y), which is formally 0(1) in the large-/?o limit. However, since 
b(mw) ~ 0.07 is numerically small, it will be a good approximation to keep only the 
one-loop contributions to the matching conditions. They read (in the NDR scheme) 

C ( 2\m w ) w y b(m w ) , Cf\m w ) w - E (x t ) b(m w ) , Cj(m f ) w F(x t ) , (13) 

where x t = {m t /mw) 2 , and 

~ /N 2, . x(18-llx-a; 2 ) x 2 (15 - Vox + Ax 2 ) , 

g oM = - 5 (1 +!■»)+ 12(1 _ j)3 ' + 6(1 _ j)4 ; ln,, 

s(2 + 5a; - x 2 ) 3a; 2 
F( ^- 8(1 i|T^ ln - < 14 > 

In Table |IJ we compare our results for the Wilson coefficients obtained at subleading 
order in the large-/?o limit with those derived using the standard RG evolution at NLO. 
Given the great simplifications that occur in the large-/?o limit, it is remarkable that the 
numerical values obtained in the two schemes are very close to each other. 
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(c) (d) 

Figure 1: Representative examples of contributions to the hard-scattering kernels 
in the QCD factorization formula including bubble resummation of the gluon prop- 
agators: (a) vertex contribution, (b) hard-scattering with the spectator quark, (c) 
penguin contribution, (d) chromo-magnetic dipole contribution. 



3 Borel resummation for the decay amplitudes 

Having obtained expressions for the Wilson coefficients in the large-/5o limit, our next 
task is to calculate the hard-scattering kernels in the QCD factorization formula using 
the same resummation scheme. There are three different types of kernels, correspond- 
ing to "non-factorizable" vertex corrections, "non-factorizable" gluon exchange with the 
spectator quark in the B meson, and penguin contributions. Representative examples of 
each kind are depicted in Figure [I], where (c) and (d) are both part of the penguin con- 
tribution. The weak decay amplitudes are obtained by evaluating the hadronic matrix 
elements of a "factorized transition operator" consisting of products of the parameters 
a>i and two quark currents, whose matrix elements are known in terms of meson decay 
constants and semi-leptonic form factors [TJ, |3|. As an example, we quote the resulting 
expressions for the B — > tttt decay amplitudes. They are 



A(B° -> 7T+7T 



A„ ai+ X P (al + 



r X a 6J 



p=u,c 



A-r, 



-V2A(B -> 7T 7T°) = \ u (ai + o a ) A nn , 



where 



V2A(B- -> 7T-7T ) - A(B° -> 7T + 7T- 



Anr = ^ « - ml) F^{ml) / w 



2m? 



m b (m u + m d ) 



(15) 
(16) 
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and we have neglected small electroweak penguin contributions and power-suppressed 
weak annihilation terms. Similar results hold for other decays such as B — > ttK. The 
general expressions for the a, parameters at leading order in the heavy-quark expansion 
read (we quote results for B — ► ttK decays; those for B — ► nir are obtained by replacing 
K — > 7r everywhere) 



oi = Ci(/i) + 

a 2 = C 2 (/i) + 

04 = C 4 (/i) + 

«e = CeO) + 



N 

CM 

N 

CM 

N 

CM) 
N 



1 + 



CfoM T/ / n , C F na 3 (n) 
— ^ — Wjf(/iJ H ^ BkM) 

C f ol s {ijl) 



Air 
C F a s (fi) 



N 

C F na s (n) 
V*W + — ^ ^(/i) 



6- 



47T 

C F a s (/j) 
Air 



C F TiaM 
VkW + 77 H K M) 



N 



+ 



47T 



AT 



47T 



N 



(17) 



The quantities V, H and P include the NLO corrections from the vertex, hard-spectator 
and penguin diagrams, respectively ||. The subscript n 
twist of the corresponding light-cone distribution amplitude of the kaon. 



2, 3 on Pk u refers to the 



3.1 General results for the Borel transforms 

The summation of terms of order /?o _1 a" in the kernels V, H and P is achieved by 
repeating the calculations of || inserting an arbitrary number of fermion loops on the 
gluon propagators (see || for details). The results are then written as integrals over 
the Borel variable u. Using the fact that a s = 0(1/ (3q) and the explicit results for the 
Wilson coefficients obtained in the previous section, it follows that 

a 1 = l + 0(1/(31) (18) 

becomes trivial in the large-/?o limit. Also, the parameters 03 and 05, which do not 
enter the B — > tttt decay amplitudes in flTSP, vanish to this order, i.e., (13, 05 = 0(1/ (3$). 
These results compare well with the full NLO analysis in [[3], where it was found that 
a\ ~ 1.03 + 0.02z, while the values of 03 and 05 are below 1% in magnitude. 

The results for the three remaining coefficients are non-trivial. In the case of a 2 , both 
the vertex term V n and the hard-spectator term H n x contribute. The vertex term has 
been computed in the large-/?o limit in H, with the result that 



4vr fa 



r *(/*) 




3 ) + 3 i n + f du e-VHm,) e 5«/3 Sv , u ) 

/ b(m h ) J 



x ' v ' 



(19) 



and the Borel transform 

1 



S v (u) = fdx$Jx) -e 5m/3 - y — — y-^- — -(-Kcotnu + in 
J \u (1 — u)(2 — u) 
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(3-u)T(2-2u)T(u) 
r(3-u) 



dz 



[z(x + (1 — x)z)Y 



(20) 



Here and below we take advantage of the symmetry of the pion distribution amplitude 
under x <-> (1 — x) to simplify the results. The term proportional to e -5 "/ 3 subtracts 
the pole at u = 0, as appropriate in the MS scheme. In the large-/?o limit, the hard- 
spectator contribution is obtained from bubble resummation on the gluon propagator of 
a tree diagram, so the result is scale independent and given by 



7ra s (/i) 
N 



H.M - I" / due""^ e 5 "/ 3 S H (u) , 
Po 



(21) 



where (we use the short-hand notation y = 1 — y) 

i i i 

S H (u) 



A 



(0) 



y c?x y y ^ 



$ ff (z) $ fl (£) 



x(yO 



l+u 



(22) 



Combining the various pieces, it is readily seen that the scale dependence cancels between 
the Wilson coefficients and the kernels. The RG-invariant result is 



a 2 



1_ 2Cp_ 
N + A% 



dfe / 7 (6) 



26 



+ / dne-"/^") e 5 "/ 3 \s v {u) + S H (u)} + O(l//3 2 ) 



(23) 



Let us now discuss the resummation for the penguin coefficients a\ 6 . It follows from 
our earlier results that in this case the vertex and hard-spectator contributions enter 
only at 0(1/(3$), and so it suffices to calculate the penguin kernels P£ n in (pf). The 
cancellation of the scale dependence between the Wilson coefficients and these kernels is 
subtle and relies in an essential way on our treatment of the Wilson coefficients in the 
large-/5o limit. Consider the coefficient a\ as an example. For fixed momentum fraction 
y, we findQ 



P a p 4 (y) = C[ 1) (fi) + 



N 



Cf 
N 



+ / du e - u/b{mb) e 5u/3 (-y- ie)' 



x 



(3o 



CPW + CPO*)] (^\n^-G(0,y) 



(o) 



G(s p ,y) __§8-L + 0(l/^), (24) 

J J 

1 For the sake of this perturbative argument it is legitimate to interchange the order of the integrations 
over u and y. 



4 m b 2 

- In 1 — 

3 /x 3 
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where s p = {m p /m b ) 2 for p = u,c, and 



G(s, y) = —4 J dz z(l — z) ln[s — z(l — z)y — ie] . 



(25) 



Noting that G(0,y) = |[~ — ln(— y — ie)], replacing nf/(3 — > —3/2, performing the 
integral over u, and using the results from Section [|, we obtain 

(3 al(y) 



2 ^C?\ri + 2 ^b(-ye-^rnl-ie) 



x (/i) ( -2 In -j- + - - \n(-y - ie) I 

+ | In ^ + ^ - J G( 8p , y)-^-) + 0(l//3b) 
3 /J 3 2 |/ J 



(26) 



2C P 



Ci 1} (M) + 6(M) 



1 1 



1 



- + ~G(0,y)-~G(s p ,y) 



(0) 

89 



+ 0(1/ p ) 



where M 2 



— y e 5//3 m 2 . — ie, and in the last step we have used the evolution equation 



for in fllT]) . It follows that 04 (y) is indeed independent of the scale /x. The terms 
proportional to Cj 1 * 1 (//) + Cg (/i) in (p4|) are absorbed by setting the scale in the Wilson 
coefficients to be the (time-like) momentum M carried by the gluon propagator. Having 
established that the resummed expression for a\ is scale independent, we return to the 
original expression in ([24]), set /i = m b , and integrate over the light-cone variable y to 
obtain the final results 



p 



2C F 
Wo 

2C F 
Wo 



C| 1} (m b ) + jdu e - u/bimb) e 5u/3 S P ^(u) 


00 

C ( l\m b ) + / due- u/h{mb) e 5u/3 S p P3 {u) 



+ O(l/f3 2 ) 



(27) 



where 



Sl 2 {u) 



dy $k(v) 



gu 



^ + ^Ci 1 \m b )G(0,y)--G(s p ,y) 



1 



(0) 



dy$ p K (y)- 



y" 



\ + lc { l\m b )G{V,y)- l -G{s p ,y)-C^ 



(28) 



Equations (|I8"D , (f23|), and (p^) are the main results of this work. They allow us to 
study the effects of higher-order /3q — corrections on the weak decay amplitudes for 
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processes such as B — > tttt and B — > 71 K. Compared with the general expressions in (|T7|), 
the terms that survive at NLO in the large-/3o limit are rather simple. However, these 
terms still capture the bulk of the non-trivial strong-interaction effects in non-leptonic B 
decays. In particular, it has already been noted in 0, [| that ai is always very close to 1, 
that hard-spectator interactions are the most important for the coefficient 02, and that 
the penguin matrix elements give the dominant NLO contributions to <x\ 6 . Before we 
evaluate our results numerically, we now simplify the expressions for the Borel transforms 
by using a model for the meson distribution amplitudes. 



3.2 Results with asymptotic distribution amplitudes 

The exact integral expressions for the Borel transforms of the hard-scattering kernels 
simplify if the asymptotic forms of the pion and kaon light-cone distribution amplitudes 
are adopted. They are 



& n (x) = 6x(l — x) 



$ K (y) = &y{i-y), *£(y) = i- 



(29) 



In addition, for the hard-spectator term we employ a simple model for the leading-twist 
distribution amplitude of the B- meson suggested in [[HJ, according to which 



MO 



1 



-mo 



6 



m B 



(30) 



The hadronic parameter Ab = O(Aqcd) is defined in terms of the first inverse moment, 
Jo(rf£/0 ^b(0 = t^b/^b- Its value at present is rather uncertain, a typical range being 
A = (350 ± 150) MeV. 

The results obtained for the Borel transforms of the vertex and hard-spectator terms 
after integration over these distribution amplitudes are 



Sv(u) 



-5u/3 



67r(cot ttu + i) 
u " (1 -u)(2-u) 2 

6r(l - 2u)Y{u 



Stt(u) 



(2 - it) r(3 - it) 

18vr 2 f B f K 



1 + 2u(l - 2u) ( V(l - 2m) - V(l - u) - - 



2r(l 



u 



(m B \ u 



N m B X B F Q B - K (0) (1 u) (2 it) V \ B 



(31) 



where ip(z) is the digamma function, defined as the logarithmic derivative of the gamma 
function. The results for the penguin contributions read 



S p P Ju) 



2 3 1 

(2-u)(3-u) + 2 C * ){mb) Gk{ °> U) ~ 2 0K{Sp > U) 



6Q 



(o) 



(l-u)(2-u) 



S 



P,3 



U 



3(1 -u) + 2 C i 1} ( m b)GK(0,u) --Gk(s p ,u) 1 



a 



(0) 

89 



u 



(32) 
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where we have defined 
i 



Gk(s,u) = Jdy<5> K (y)y- u G(s,y), g K (s,u) = J ' dy & K (y) y~ u G(s,y) . (33) 
o o 

For p = u, these functions can be evaluated in closed form, yielding 

G K (0, U) = — rj- r - + «7T + + 



(2-u)(3-u) V3 2-u 3-u 

^ '") = 3(i^)(5 + iT + T^)' (34 > 

For the case of the charm-quark loop, the functions Gk(s c ,u) and Qk(s c ,u) with s c ^ 
do not have a simple analytic form. We choose instead to write them as infinite sums 
of poles in the Borel plane, by expanding the logarithm in (|25| ) in powers of y and then 
performing the integrals. We find that 

_ -41ns c ^ 0r ,_ w r(n + 2) 6 

^ l ' c ' MJ (2-«)(3-u) + ^i 2n 1 cJ r(n + f) (2 + n - tz)(3 + n - tz) ' 

ri \ -21ng c , ^ V^ (A ^ n T(n + 2) 1 

The sums converge as long as s c > 1/4 and can be expressed in terms of hypergeometric 
functions. These functions can then be analytically continued to the physical region 
s c < 1/4, where they develop a non-zero imaginary part (recall that s c — > s c — it). 



3.3 IR renormalon ambiguities and power corrections 

Expansion of the Borel transforms (times e 5 "/ 3 ) in powers of u generates the perturbative 
series of the kernels in the MS scheme. The Borel sums of these series are obtained by 
performing the M-integrals exactly. However, the Borel transforms contain poles along 
the integration contour, which render the results of Borel summation ambiguous. This 
is the well-known problem of IR renormalons, which reflects the fact that the original 
perturbation series were asymptotic. 

For a given Borel transform, the residue of the nearest pole provides a measure of the 
renormalon ambiguity, which is conventionally taken as an estimate of the power correc- 
tions that should be added to the perturbation series in order to get a well-defined result. 
We will now work out the structure of these leading singularities for the parameters aj. 
Note that the product 

~u/b(m b ) 5n/3 _ ( ^^MS 



Av\ 



2 u 



(36) 
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is RG invariant. The constant 5/6 arises in the calculation of the fermion-loop contribu- 
tion to the gluon self-energy in the MS scheme and compensates the scheme dependence 
of Ajyfg. We will thus express our results in terms of the invariant scale parameter Ay. 

If the light-cone distribution amplitudes are given by the model used in the ex- 
plicit calculations of Section |3.2| , then all of the Borel transforms have a leading pole at 
u = 1. For the hard-spectator term this leads to a power correction of order Ay/m b , 
which is expected, because the effective scale in the hard-spectator terms is not ml but 
n\ ~ Aqcd^6 U- In all other cases the leading power correction is of order (Ay/m b ) 2 . 
Expanding the Borel transforms about the positions of the poles, and including terms 
up to second order in the heavy-quark expansion, we find 



Aa 2 
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cf(m b ) 



(37) 



where A^ = (m B /m b )(Ay / X B ) = 0(Aq C d) in the hard-spectator contribution to Aa 2 . 
These expressions may serve as a model for the leading power corrections to the decay 
amplitudes arising from "non-factorizable" soft gluon exchange. Because most of the 
corrections are suppressed by two powers of A v /m b , their numerical effects turn out to 
be rather small. Only a 2 receives a first-order correction that can be sizeable. However, 
the effect of that correction can be absorbed into a redefinition of the hadronic parameter 
X B , which in any case is poorly determined theoretically. Therefore, power corrections 
from soft gluon exchange have a minor impact on the phenomenological analysis of 
charmless hadronic B decays in the framework of QCD factorization. 

Let us briefly discuss the changes that occur if the light-cone distribution amplitudes 
are allowed to have an arbitrary endpoint behavior, such as 



1 + 0(x) 



(38) 



as x — > 1 (and similarly for x — * 0). The model used above corresponds to 5^, 5k, Sb — 1 
and S K = 0. These results would continue to hold at any fixed order in a Gegenbauer 
expansion of the distribution amplitudes, but may be altered if an infinite number of 
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Gegenbauer moments is included. The effect of a more general endpoint behavior on the 
position of the poles in the Borel plane depends on the function being considered. As 
seen from the Borel transform corresponding to the vertex function Vm has a pole at 
u — 1 regardless of the endpoint behavior of the distribution amplitude of the meson M. 
For the penguin and hard-spectator terms, on the other hand, changing the distribution 
amplitudes will change the position of the leading IR renormalon poles. For the penguin 
terms the modifications are simple: the leading poles in the twist-2 and twist-3 penguin 
contributions are shifted from u = 1 to u = 5 k (twist-2) and u = 1 + 5 K (twist-3), 
respectively. The effects on the hard-scattering term are slightly more involved. If the 
5-meson and kaon distribution amplitudes have the same endpoint behavior, then there 
are both single and double poles at u = 5 k = 5b- On the other hand, if 5b differs from 
5k, then the position of the leading simple pole is determined by whichever is smaller. 



3.4 Fermion loops beyond the large-/3 approximation 

Whereas in most applications of the large-/5o limit the replacement rif — > —3/3q/2 after 
the calculation of fermion bubble graphs can be motivated by a diagrammatic analysis, 
this replacement is hard to justify for the penguin coefficients entering the analysis of 
non-leptonic B decays. The fermion-loop contributions proportional to rif are usually 
accompanied by gauge boson and ghost loops. Calculations using bubble resummation 
are often good approximations to full QCD because the replacement rif — > —3(3q/2 
(approximately) accounts for the extra diagrams not related to fermion bubbles. The 
factor of nf associated with the penguin coefficients is not of this nature. It arises 
from the contraction of two fermion fields of a 4-quark operator in the effective weak 
Hamiltonian. Such fermion loops are not accompanied in an obvious way by other QCD 
loop diagrams.^ To produce reasonable numerical results for the Borel-resummed kernels 
will require moving beyond the large-/?o approximation. 

A proper counting of fermionic penguin contractions can be achieved a posteriori by 
making the replacement 

I cf(m & ) G(0, y) - W n _ 2 ) G(0, y) + G(s c , y) + G(l, y)} (39) 

in the expressions for the Borel transforms S P2 and S P3 in (p8|). Upon setting rif = 5, 
(5^) accounts for loop contributions from three massless quarks, the charm quark, and the 
bottom quark (cf. eqs. (49) and (54) in 0). Formally, the two sides of (|39D coincide in the 
large-/?o limit, and so making this replacement is equivalent to adding some higher-order 
terms to the renormalon calculus. Since there are other terms of subleading order in the 
large-/?o limit not taken into account, this procedure may seem arbitrary. Nevertheless, 
the above replacement is physically motivated and, in particular, reproduces the proper 

2 Such diagrams would appear at two- loop order when one calculates matrix elements of the chromo- 
magnetic dipole operator. We have performed a partial analysis of such higher-order contributions and 
found their effects to be strongly suppressed. 
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imaginary parts and heavy-quark mass dependence of the fermion loops. From now on, 
we will refer to the scheme defined by ( p9[) as the "modified" large-/?o limit. 

A replacement analogous to fl55|) must be made in fl32|). The resulting expressions 
for the leading IR renormalon ambiguities in the coefficients Og also change. In the last 
two equations in (|37|), we must replace 

(n/ _ 2) ( 21n -^) + 2 ^]. (40) 
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4 Numerical results 

In addition to providing a model for power corrections, the results of the previous sections 
allow us to study the numerical significance of higher-order perturbative contributions 
to the hard-scattering kernels in the large-/?o limit, which may be considered as a toy 
model for the asymptotic behavior of perturbation theory. Partial perturbation series 
for the quantities dj, including terms up to order Py^a™, are obtained by retaining only 
terms up to 0(w n_1 ) in the Taylor expansions for the Borel transforms about u — 0. The 
resummed series (in the large-/3o limit) are obtained by performing the Borel integrals 
exactly, regularizing the singularities along the integration contour using the principal 
value prescription. The sum of the residues of the poles at half-integer values of u provides 
an estimate of the irreducible ambiguity in the definition of the resummed perturbation 
series. Whereas only the leading poles were considered in the previous section, exact 
expressions for the sum of all residues are used in our numerical analysis. 

In Table ||, we show different perturbative approximations to the coefficients a<i and 
(a p A + r x ag) entering the B — > mr, %K decay amplitudes. The first row gives the results 
obtained using the standard QCD factorization approach described in ||. The following 
rows show the results obtained in the large- /3q limit, at one-loop and two-loop order, and 
from a numerical evaluation of the Borel integrals. The last row contains the renormalon 
ambiguity. For the penguin coefficients, the default entries refer to the modified large-/?o 
limit, whereas results obtained in the strict large-/5o limit are given in square brackets. 
The only sizeable renormalon ambiguity arises from the hard-spectator term, which 
contributes to the real part of a 2 and corresponds to a first-order power correction from 
soft gluon exchange with the spectator quark. The ambiguity in the penguin coefficients 
is of order (Aqcd/^?)) 2 and numerically very small. Based on our renormalon analysis, 
we thus find no evidence for a large dynamical enhancement of penguin amplitudes, in 
contrast with some recent conjectures |20|, j2~T| . 

In regards to the numerical significance of higher-order perturbative corrections, the 
most striking effect is the large increase of the imaginary part of the coefficient a 2 as one 
goes to higher orders in perturbation theory. An analysis using momentum distribution 
functions shows that the imaginary part of the hard-scattering kernel contributing to 
a 2 is dominated by contributions from scales of order a GeV, indicating that it can 
receive large higher-order corrections M. The corresponding strong-interaction phase 
may therefore be larger than that obtained from the NLO analysis. 
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Table 2: Numerical results for the combinations of dj parameters entering the 
B — ► 7T7T decay amplitudes. Input parameters are f n = 131 MeV, fs = 180 MeV, 
F(f^(0) = 0.28, \ B = 350 MeV, m b = 4.2 GeV, m c = 1.3 GeV, and r x = 1.197. 
Perturbative results refer to the scale /i = m b in the MS scheme except for the hard- 
spectator contribution to a 2 , where the scale is \i = v^W^fe with = 0.5 GeV 0. 
Borel resummed results are scale independent. Values given in brackets refer to the 
strict large-/? limit. 





a 2 


a\ + r x al 




NLO 


0.086 - 0.084z 


-0.081 -0.034z 


-0.091 - 0.0142 


1-loop 


0.084 - 0.0752 


-0.082 - O.O262' 
[-0.121 -0.071i] 


-0.091 - 0.007i 
[-0.130 - 0.0532] 


2-loop 


0.099 - 0.112i 


-0.072 - 0.0362' 
[-0.107 -O.lHi] 


-0.090 - 0.0132 
[-0.125 - 0.0912] 


Borel sum 
Ambiguity 


0.094- 0.1682 
0.058 - O.OO82' 


-0.062 -0.0312 
[-0.070 - 0.115i] 

0.0004 [0.002] 


-0.088 - 0.012i 
[-0.096 - 0.096i] 

-0.0004 [0.001] 



In the strict large-/?o limit (shown by the numbers in brackets), the penguin coeffi- 
cients a\ and ag mirror the behavior of a<i- On the other hand, when the physical value 
nf = 5 is used as in (|3"9"|), the phases are reasonably close to those at NLO and do not 
receive large higher-order corrections. This shows that after the replacement in (j39|) 
there is a strong cancellation between the different terms in (28), which is maintained 
at higher orders. While the replacement in ( j39|) was motivated on physical grounds, it 
is an open and important question whether such a delicate cancellation will persist once 
the full set of NNLO corrections is included in the theory of QCD factorization. 



5 Phenomenological applications 

This section addresses the effects of renormalon ambiguities and higher-order corrections 
on the phenomenology of B decays into two light pseudoscalar mesons. The results 
of Section |] have shown that the most notable effects are the ambiguity in the hard- 
spectator term and an increase in the imaginary parts of the kernels at higher orders in 
perturbation theory. We limit our phenomenological analysis to quantities that highlight 
these effects in a transparent way. 

The B ± — > tt^tt decay amplitudes depend strongly on a 2 (see (|i~5"D) and are therefore 
a good probe of the renormalon ambiguity in the hard-spectator term. Using the data 
in Table |2], and taking [Vub/Kfel = 0.085 and |V^| = 0.041, we find for the CP-averaged 
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Figure 2: Predictions for the B — > 7r ± i ; r T and B T — > 7r i^ T CP asymmetries in the 
strict (right) and the modified (left) large-/?o limit. The different curves refer to the 
results at NLO (solid), one- loop order (long dashed), two- loop order (short dashed), 
and Borel-resummed (gray bands). For the Borel-resummed results, the width of 
the bands shows the renormalon ambiguity. 

branching ratio 

10 6 Bt(B ± -> 7T ± 7r°) = (4.57; 4.73; 4.74 ± 0.50) , (41) 

where the numbers in parentheses refer to one-loop order, two-loop order, and Borel- 
resummation, respectively. The NLO result obtained using QCD factorization is 4.80. 
(The central value quoted in || is about 10% larger due to the inclusion of electroweak 
penguin contributions.) The magnitude of the renormalon ambiguity may be compared 
with the error estimate given in [Q, where an uncertainty of i^l (pars.) ±0.3 (power) was 
assigned to account for the effects of parameter variations and power corrections from 
higher-twist distribution amplitudes. Since this is the mode that is the most vulnerable 
to the ambiguity in <Z2, we are led to the conclusion that no large power corrections from 
soft gluon exchange are present in B — ► tttt, ttK decays. 

The CP asymmetries of certain decays can be very sensitive to the imaginary parts 
of the hard-scattering kernels and show clearly the effects of Borel resummation. As 
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examples, we discuss the asymmetries for B — > n^K^ and B — > tt°K t , which to first 
approximation are given by 

A CP (B° -> 7r + K~) oc - sin 7 • lm(a c 4 + r x a c 6 ) , 

A CP (fi- - tt ^-) oc -sin 7 • Im f "^^ ) , (42) 

where ~ F^ K (0) / (0) x f w / fx ~ 0.9. Using asymptotic distribution am- 
plitudes, the ctj parameters for B — ► 7T.K" decays are identical to those for I? — > 7r7r. 
In Figure |2], results for the two CP asymmetries are shown both in the strict large-/9o 
limit (right) and in the modified scheme discussed in Section |3.4 In both schemes, the 
asymmetries increase strongly at higher orders in perturbation theory, such that the 
Borel-resummed results are almost a factor of two larger than the results found at one- 
loop order. The two-loop contributions appear to be the most important higher-order 
effects. This trend suggests that the asymmetries predicted from exact NLO calcula- 
tions in || will become larger upon the inclusion of higher-order terms. (As concerns 
the magnitude of the CP asymmetries, we believe that the true values will be some- 
where in between the two schemes.) It is also apparent from the figure that, in all cases, 
higher-order perturbative corrections are much more important than the power correc- 
tions suggested by the values of the renormalon ambiguities. This provides additional 
motivation for calculating the imaginary part of the hard-scattering kernels at two-loop 
order. In fact, the two-loop imaginary parts (combined with the one- loop real parts) of 
the decay amplitudes are required for a consistent description of the CP asymmetries at 
NLO in RG-improved perturbation theory. 

6 Conclusions 

In this paper we have used the renormalon calculus and the large-/?o limit of QCD to 
develop a toy model for higher-order corrections to the QCD factorization approach as 
applied to B decays into two light pseudoscalar mesons. We have presented novel results 
for the Wilson coefficients in the large-/3o limit, and have used them to show that the 
RG invariance of the hard-scattering kernels entering the QCD factorization formula is 
preserved under Borel resummation. Both the Wilson coefficients and the coefficient 
functions Oj appearing in the amplitudes for B — > im, ttK decays simplify greatly in the 
large-/3o limit, but still allow for a realistic study of power corrections and higher-order 
perturbative contributions to QCD factorization. We have presented exact expressions 
for the Borel transforms of the hard-scattering kernels in the factorization formula, from 
which is it possible to derive the terms of order Py^a™ in the perturbative expansions 
of the kernels. We have analysed the power corrections from soft gluon exchange and 
obtained a simple model for the first and second-order corrections to the parameters 
dj. The only numerically significant power correction arises from soft gluon exchange 
with the spectator quark, which introduces an error in the CP-averaged 5 ± — > 7r ± 7r 
branching ratio of about ten percent. This is roughly the same as the errors associated 



17 



with the hadronic uncertainties and power corrections that were estimated in [^]. The 
power corrections to the real part of the penguin coefficients are of order (AQ CD /m;,) 2 , 
and those to the imaginary part are zero within the large-/5o limit, which offers numerical 
evidence against models proposing dynamically enhanced penguins. 

A recurring result of our study is that the imaginary parts of the hard-scattering 
kernels may be larger than expected from previous applications of QCD factorization. 
A calculation of the imaginary part of the kernels at two-loop order would be highly 
desirable, since it would give an answer to the question of whether there are large multi- 
loop contributions to these decays that are not captured in the large-/3o limit. 
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